Generalized grand-canonical ensemble theory for interacting Bose-Einstein systems 
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We use the maximum information principle to include particle-interaction into the grand-canonical 
theory of BECs. The inclusion of the particle-interaction elucidates why thermodynamic calculations 
for BECs by the grand-canonical ensemble for the non-interacting case are in coincidence with up to 
date experimental results. However, in our generalized theory we can show that a BEC experiences 
a real and abrupt phase-transition in contrast to the smooth phase-transition predicted by the 
non-interacting grand-canonical ensemble. In addition, we discuss possible effects due to a thermal 
environment and the experimental probing from a general point of view. 



PACS numbers: 34.50.Dy, 03.65.-w, 03. 75. Be 



I. INTRODUCTION 

The prediction of a macroscopical, condensed state for 
bosonic systems in P, [3] started an extremely fruitfuU 
theoretical development in order to describe this phe- 
nomenon in all details. The theory of Bose-Einstein con- 
densation has been worked out macroscopically within 
the framework of the canonical and grand-canonical en- 
semble, e.g. [H, 0, H, 0, but only for cases without 
interaction between particles. Particle-interaction due to 
s-wave scattering, on the other hand, was included in 
microscopical approaches related to quantum field the- 
ory at zero temperature. Amongst the microscopical ap- 
proaches, the standard references are the theory by Bo- 
goliubov [il for the quasi-particle picture and the work of 
Gross fi\ and Pitaevskii [sj for the inclusion of particle- 
interaction into the framework of the Schrodinger theory. 
Contact between the interacting field-theoretical and the 
non-interacting ensemble approaches is usually made by 
including microscopically calculated spectral corrections 
into the statistical description of the many-particle sys- 
tem by various techniques that are reviewed in, e.g.. 

Our work will pursue a completely different, exclu- 
sively macroscopical ansatz for an equilibrium ensemble, 
that also may be called a holistic one. Our ansatz is based 
on the well defined fact that all available and necessary 
information about a macroscopical physical system can 
be obtained from its entropy. The meaning of entropy as 
a measure for information was first recognized by Shan- 
non an d was brought back to its physical context by 
Jaynes |l2j | . In a variational principle was developed 
that allows us to derive the relevant physical distribu- 
tion functions under several constrains imposed on the 
entropy by the method of Lagrange-parameters. This 
variational principle is also known as the maximum in- 
formation principle (MIP), e.g. [l^. By the MIP, all 
statistical ensembles known are derived most easily when 
the appropriate constrains are applied. 

In our approach we will use the MIP to include 
the particle-interaction for the BEC-case into the sta- 
tistical grand-canonical ensemble description in a non- 



perturbative way. The generality of the MIP ensures, 
that only the most elementary knowledge about the 
BEC-system is required: the Hamiltonian, the particle- 
conservation and the symmetry of the many-particle 
wave- function. 

We further will discuss environmental effects onto the 
generalized BEC-ensemble, that are directly related to 
its non-linear structure. We will treat the environment 
in the sense of equilibrium statistics which is valid for a 
system after cooling-down and relaxation to equilibrium 
have come to an end and only quasi-stationary interac- 
tions between system and environment are left. 



II. BOSE-EINSTEIN-STATISTICS INCLUDING 
PARTICLE-INTERACTION 

The quantum system describing interacting bosons in 
terms of the density-matrix is given by the Hamiltonian 



H = tr[H] = tr[HoP 
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where Hq = T + U is the non-interacting part of the 
system containing the kinetic energy T and the external 
potential U. 

We abbreviate the entropy by S* = tr['S'] = 
— tr [/51n[/5]] and together with Eq. ^ we obtain the 
variational equation for the density-matrix of the grand- 
canonical ensemble, 
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The Lagrange-parameters Xi in Eq. ^ ensure that the 
system is constrained with respect to the conservation of 
the energy and the particle number. Variation of Eq. Q 
with respect to p yields 



\n[p] + Ai [Ho + 2gp) + X2 N ^ 



(3) 



Equation ([3]) is solved in terms of the Lambert function 
A, which is related to the product-logarithm. In terms of 
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the Lambert function A the density-matrix of the inter- 
acting system is readily given by 



= A(/5o;/3<?) = 



i^Pgy-'pi 
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From the laws of thermodynamics we find the Lagrange- 
parameters with their familiar meaning Xi = P — 
l/kgT and A2 — (3 ^i, leading to the well-known density- 
matrix of the non- interacting system po = exp[— /? (Hq — 
fJ-N)]. The symbolic equation Q becomes meaningful 
when we calculate the partition function by applying the 
trace-operation 
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In Eq. ([5]) we have defined partial partitions 



/3(fp-A') 



(5) 
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The partial partitions in Eq. ^ describe the degree of 
the interaction of the system according to the power-law 
expansion of the Lambert function Eq. ([4]). Quite sur- 
prisingly, the partial partitions have the same structure 
as in the non-interacting case, revealing the presence of 
particle- interaction only through the additional factor j. 
This special property is a direct consequence of the con- 
tact interaction assumed to be valid in the BEC case. 
When interactions are non-local the density-matrix is 
also non-local, inducing tremendous complications to Eq. 
([4]) . The special property of Eq. ([6]) can be understood 
by recalling the general property of the density-matrix, 
namely 



tv[p] > tr[p2], tr[p]2 > tr[p2] 



(7) 



The inequalities Eq. ([7]) follow directly from the Cauchy- 
Schwarz inequality and can be extended to general pow- 
ers. Thus, the partial partitions Eq. © obey the relation 
Zj > Zj+i. 

We rewrite the partial partitions Zj by introducing 
partial grand-canonical potentials $j 



3\ ' 



Zj = cxp[-/3$ 
p 



(8) 



which will prove useful for the following calculations. We 
proceed by calculating the grand-canonical potential S = 
— /3^^ In [-2], reading 



-1 ^ (-If 



k+l 



with 



k=l 



E^^T-^(2/35r^exp[/3A<I>iJ I ,(9) 



where we have separated the non-interacting part giving 
the first term and have already expanded the remaining 
logarithm into a power-series with the convenient nota- 
tion A$ij = - ^j. 

So far, our calculations have been carried out on a 
level of full generality. This allows us to apply our above 
theory to any bosonic quantum system where the free 
spectrum ep or, equivalently, the free density of states 
D{e) is known. 

A system of considerable interest in the context of 
BECs is the harmonic oscillator. In three dimensions, 
its density of states is given by [1| , 



Die) 



(10) 



where in Eq. (jlOp Q is the geometric mean of all frequen- 
cies u) — (wi wa)^/"^. When Eq. pH)) is inserted in 
Eq. ([8]) and the sum in Eq. ([8]) is taken in its integral 
limit which applies for large particle numbers, the partial 
grand-canonical potential is easily calculated, giving 



(11) 



where the symbol Li,i in Eq. (|lip is the usual notation 
for the poly-logarithmic function of order n. Knowing 
the partial partitions explicitly, the occupation number 
of the excited states can be obtained by the laws of ther- 
modynamics 
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In Eq. C(3) is the Riemann zeta function, and the 
function x describes the effect of the particle-interaction 
on the occupation of the excited states. In Fig. 1 we have 
plotted the behaviour of Eq. (fT^ and its non-interacting 
counterpart in arbitrary units for comparison. Clearly 
visible is a sudden fall-off of the number of occupied ex- 
cited states at a certain temperature, that we will call 
the condensation temperature Tcond from now on. When 
the system has reached this temperature from above, all 
excited states experience a sudden deoccupation, such 
that the only occupied state left is indeed the ground- 
state of the system. This special behaviour of Eq. (fT^ 
reveals the Bose- Einstein phase-transition for interacting 
systems not to be smooth. Our result is in sharp con- 
trast to the Bose-Einstein phase-transition predicted for 
non-interacting systems, where the occupation number of 
excited states only vanishes as the temperature T reaches 
zero, e.g. |3„ Strictly speaking, such behaviour does 
not contain any kind of phase-transition. 

From the behaviour of Eq. in Fig. 1 it becomes 
clear that, when Tcond is known, the occupation number 
is sufficiently described by a step-function 



W=2 



C(3) e[T -T, 
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FIG. 1: Occupation numbers A^cx for all excited states for the 
interacting case (full line) and the non-interacting case (dashed 
line). To illustrate the influence of the particle-interaction, we set 
kB = h = Hi = 1 and choose g = 0.3. The temperature T^ond 
marks the point where the BEG phase-transition occurs. 



FIG. 2: Condensed fraction of particles for the interacting case 
(full line) and the non-interacting case (dashed line) in arbitrary 
units for g = 0.3 and a critical temperature Tc = 5. The vertical 
line marks the condensation temperature T^ond- 



It is clear from Eq. p3)l that the existence of the con- 
densation temperature Tcond does not harm the meaning 
of the critical temperature T^, see e.g. [H, as it should. 
The critical temperature remains to be the temperature 
at which, when approached from above, the ground-state 
of the system starts to get occupied. As a matter of 
fact, this occupation of the ground-state does not depend 
on the presence of particle-interaction and is a pure ef- 
fect of cooling down a bosonic system. This explains 
the success of the non-interacting grand-canonical en- 
semble in this point. Only when Tcond is reached the 
particle-interaction gives rise to a sudden completion of 
the ground-state occupation. Thus, one could speak of a 
sliding phase-transition in the BEC-case, starting at the 
critical temperature Tc and finishing at the condensation 
temperature Tcond- 

According to Eq. the law for the condensed par- 
ticle fraction stays in a form similar to the known one. 



No T3 



(14) 



The influence of Tcond can be seen from Fig. 2, where we 
have plotted the condensed particle fraction illustratively 
in arbitrary units. For temperatures Tcond < T < Tc 
the behaviour of Eq. shows the well-known power- 
law tail, e.g. [5|], but for temperatures < T < Tcond, 
according to Eqs. p2l[T3)l the number of excited particles 
falls off to zero and every particle is to be found in the 
ground state. 

We now discuss our results in their relation to exper- 
iments that realize BECs. The influence of the particle- 
interaction in Eq. is described by the ratio g/{kB T). 



When we scale the system in units of its critical tem- 
perature Tc by introducing T = tTc we can define an 
effective coupling constant 



9cs 
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kn Tc 



(15) 



When Eq. p5|) is inserted into Eq. p2)) . we can cast it 
into 



C(3)x(r^r3;M), 



(16) 



With the value of gcff from prominent experiments, e.g. 
[ii, m, H [S [ll, we find the interesting fact that all 
experiments have in common an extremely small value 
of gcs- For a quantitative illustration we take the JILA 
experiment 16], where ^"^Rb atoms where Bose-Einstein 
condensed. The critical temperature in (l6j was deter- 
mined to be Tc = 280 nk and the s-wave scattering 
length is found in [Tol, aint = (103 ± 5) x OBohr- From 
this data follows a critical thermal energy fee Tc = 3.86 x 
10^'^'^ J and an interaction energy oi g = 5.27 x 10^^^ J. 
The effective coupling thus gives gcS = 1.36 x 10"^-^, 
which is extremely small. Only when the dimensionless 
temperature r reaches a magnitude of r « 10""'^^ — 
10~^^, the effect of particle-interaction on the occupa- 
tion number becomes visible. Thus, in up to date exper- 
imental setups, the thermal energy at the critical tem- 
perature Tc exceeds the interaction-energy g by orders 
of magnitude. This makes the condensation temperature 
Tcond to lie only slightly above the absolute zero-point 
and the gas behaves effectively as a free gas with negli- 
gible influence of its internal interaction as long as Tcond 
close to the absolute zero-point has not been reached. 
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This remarkable fact is the reason why, besides tiny 
corrections due to finite size effects, e.g. [l^, no seri- 
ous deviations - clearly distinct from systematic errors 
due to the measurement technique - from the predic- 
tions of the free grand-canonical ensemble have been de- 
tected in experiments so far. From technical remarks on 
the experimental feasibility, e.g. p^ . concerning espe- 
cially measurements at very low temperatures close to 
the absolute zero-point we conclude, that the strongly 
increasing signal-to-noise ratio in this regime is respon- 
sible for the experimental inaccessibility of the extreme 
low-temperature region. However, a measurement of the 
condensed fraction of particles in a system that is close 
to a Feshbach-resonance, increasing g by orders of mag- 
nitude, should make the BEC phase-transition become 
visible as predicted by Eqs. (fT^[H|) . 



III. INFLUENCE OF THE ENVIRONMENT 

Above we have mentioned experimental limitations for 
the access of the regime of extreme low temperatures. 
As no macroscopical system can be regarded as isolated, 
we proceed now by including effects, that may arise by 
the coupling of the condensate to a thermal environment, 
where at least black-body radiation is present. It is com- 
mon sense in the BEC-field, and confirmed by experiment 
so far, that black-body radiation docs not affect the con- 
densation process, e.g. [l0|. However, this changes when 
effects of non-linear ensembles are to be studied, as we 
will show now. Our ansatz follows the spirit of typical 
system-bath models, e.g. [HI, [l^, [2^ , but will be adopted 
to the requirements of our system. The simplest possible 
interaction between the system (1) and the bath (2) is 
given by a Hamiltonian 



int = K Pi /02 



(17) 



The coupling-constant k has the meaning of some 
transition-rate. The full Hamiltonian is given by 



H 



Ho + 9 Pi] + H2 + npi p2 



(18) 



When we discard the back-reaction from the system on 
the environment, we obtain the effective Hamiltonian by 
tracing out all contributions of the bath, giving 



Ho + E2 



r(2) 



V 



z. 



(1) 



Pi -(19) 



The thermal energy and the partition-functions of the 
black-body radiation can be found in any textbook on 
statistics, e.g. We have i?2(T2) = 4:Cr/cT2, where 



a is the Stefan-Boltzmann constant and c is the speed of 
light. Further, from our calculation above, Eqs. ([51 [TT|) . 



we can get for the partitions 



Z2 — exp 



-7(2) 

Z2 — exp 
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From the partitions Eq. (PO)) it can be seen that for 
T2 ^ the influence of the environment on the coupling 
vanishes. Finally, we introduce an effective coupling 



r(ff,T'2,K) = g + R{T2)k , 
where we have defined the ratio 



RiT2 



r(2) 



Z. 



(1) 



exp 
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(22) 



From Eq. ([22)1 it can easily be read off that the coupling 
is only influenced by the environment for small, but not 
too small temperatures T2. For higher temperatures Eq. 
([^^ declines fast and renders the system with the bare 
particle-interaction only. With Eq. ([^Tjl the density- 
matrix including environmental effects, according to Eq. 
([4]), readily follows 



i-j) 



J'- 



Pi 



(23) 



where the density-matrix of the free system now shows a 
thermal energy-shift 



Po 



exp 



Pi{Ho + E2{T2) - piNi 



(24) 



It is worth noting, that this thermal shift is present even 
without any explicit coupling between system and bath. 
According to Eq. ^ we find the partition-function 

^1 = ^^^(2 Pi ry-' Z, exp[-j /?i E2{T2)] . 



J'- 



(25) 

In our following discussion, wc will focus on the influence 
of the thermal energy-shift, since the coupling has only 
a marginal effect. From Eq. (|25p it can be immediately 
read off that the thermal energy-shift, for high environ- 
mental temperatures, beats down the interacting part of 
the partition-function such that only the non-interacting 
part for Zj = 1 survives. The destruction of the non- linear 
behaviour is initiated when the exponent in Eq. ()25[) be- 
comes equal or larger than unity. Thus we obtain the 
condition 



To 



4a 



Ti 



(26) 



When Ti is taken as the temperature of the system 
Ti ^ 10~^ K, the highest allowed value for T2 is of the 
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order 10~^ K. From this rather strict condition follows, 
that the emergence of a non-linear ensemble does not 
only require a tuned bare interaction g, but also a cold 
environment. Tuning the bare interaction by a Feshbach- 
resonance without an appropriately prepared environ- 
ment, one would gain nothing. The occupation-number 
of the system above Tcond is only weakly influenced by 
the thermal energy-shift, making an interacting system 
to behave like a non-interacting system. The result of 
this is well-known: condensation is completed only when 
absolute zero temperature is reached. Thus, the thermal 
energy-shift destroys the existence of the condensation 
temperature Tcond and shifts it to zero. 

The subtle influence of the black-body radiation on 
the structure of the ensemble is only one effect of a dirty 
vacuum. Other dirty vacuum effects are discussed in [10| 
as miscellaneous small effects and all of them are under 
suspect to simply destroy the non-linear structure of the 
BEC-ensemble and let it appear as being the well-known 
linear grand-canonical ensemble. 

IV. CONCLUSION 

We have used the maximum information principle 
(MIP) to include the contact-interaction present in a 
BEC-system into the statistical description of the grand- 
canonical ensemble in a non-perturbative way. The MIP 
provides a very general approach, because it requires 
only the most elementary knowledge of the properties 
of the system: the Hamiltonian, the condition of the 
particle-conservation and the symmetry of the many- 
particle wave-function. Our analysis of the occupa- 
tion number revealed the existence of a Bose-Einstein 
phase-transition occurring at a condensation tempera- 
ture Tcond- At Tcond > all excitcd levels become de- 
occupied and the occupation of the ground-state is com- 
pleted. The existence of the condensation temperature 
Tcond does not affect the meaning of the well-known criti- 
cal temperature Tc, where, when approached from above, 
the occupation of the ground-state is initiated. The ex- 
istence of both temperatures states a sliding character 



of the Bose-Einstein phase-transition, starting at Tc and 
being abruptly completed at Tcond > 0. 

We have discussed that the experimental state-of-the- 
art is unable to access the extreme low-temperature 
regime where deviations from the free grand-canonical 
ensemble description become visible. The effective cou- 
pling constant gcs gives two possible paths of how the 
extreme low-temperature regime may become accessible 
in experiments. The first path may be a further reduc- 
tion of the critical temperature of the system, however, 
this goes necessarily along with a drastic reduction of the 
total number of particles, that can only be avoided by an 
extremely flat trapping potential. The second path may 
be a tuning of the system by a Feshbach-resonance to 
manipulate g instead of Tc. 

We also have discussed possible effects of the environ- 
ment in the case of a dirty vacuum, that is poisoned by 
black-bod y ra diation or other miscellaneous effects as dis- 
cussed in [lO|. Our discussion revealed, that tuning the 
particle-interaction by a Feshbach resonance might not 
be enough, in the experiment one has also to take care 
for an almost clean vacuum around the condensed atomic 
cloud. At least for blocking black-body radiation as good 
as possible, one could think of a second atomic system 
surrounding the system of interest in order to create a 
screening effect. When such an artiflcial environment 
could be cooled down along with the condensate, thereby 
following the above relation, Eq. a realistic chance 

would emerge for experimentally probing the non-linear 
structure of the BEC-ensemble above zero temperature. 

Our result hopefully inspires future technical develop- 
ments towards the experimental accessibility of the low- 
temperature regime, because the technical handling of a 
condensation temperature Tcond above the absolute zero- 
point would allow to study condensed atomic systems 
behaving effectively as being at the absolute zero-point 
although their true temperature is still flnite. We close 
our discussion by the general remark that our analysis 
in the BEC-case could provide some similarities to the 
phenomenon of high-temperature superconductivity. Its 
realization, however, depends on future experimental de- 
velopments. 
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